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ABSTRACT 



The structure of barotropically unstable disturbances 

In the tropics is studied with a two-level quasi-geostrophic 

model. An Ekman layer is attached to the lower boundary. 

The equations are linearized and the most unstable mode is 

found numerically by use of the initial value technique. 

CISK heating, horizontal friction and shear are introduced 

to give more realistic solutions. Computations are made for 

two shear zone wind profiles. The first is represented by 

CT, = U_ = -U tanh y/y and the second is a mirror reflection 
l i o o 

of the first and is given by = -U^ . The model is run 
for 1^2 days after which growth rates are determined as a 
function of the non-dimensional wave number ky Q . The 
associated phase angles and amplitudes for the wave numbers 
with maximum growth rates are given. The solutions with 
heating do not have a strongly preferred mode of maximum 
growth. 
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TABLE OP SYMBOLS AND ABBREVIATIONS 



A g Vertical Eddy Viscosity 

A h Horizontal Eddy Viscosity 

$ Derivative of Coriolis parameter at y = 0 

Cp Specific heat at constant pressure 

f Q Coriolis parameter at y = 0 

g Gravity 

w dp/dt 

^ gz/f Q 

R Gas constant for dry air 

p Density 

T Average temperature from the standard atmosphere 

a (R 2 T/p 2 g) (8T/8z + g/c ) 

C V 2 4; 

z Height 

Q Heating added per unit mass 

u 2 2 

V Bi-harmonic operator (V ) 

a Wind inflow' angle 

C 
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I. INTRODUCTION 



Synoptic scale disturbances in the tropical atmosphere 
have been recognized for a long time (Palmer 1951 , Riehl 
195*1). This paper is concerned with wave disturbances that 
occur in the lower troposphere and makes no distinction 
between easterly waves and equatorial waves. Starr and 
Wallace (196*1) showed that in the lower troposphere these 
disturbances have a "cold core" structure; the rising air is 
cooler than the air which is subsiding. Warm core structures 
at higher levels with significant correlations between 
temperature and vertical velocity have been found by Chang 
et al. ( 1970 ) and Nitta ( 1970 ). Vertical velocities associ- 
ated with these disturbances in general are rising air to 
the east of the wave axis and subsidence to the west (Yanai 
and Nitta, 1967 ) • 

The energy source for these waves is not known. 

Baroclinic instability is unlikely because of the smallness 
of the Coriolis parameter. Charney ( 1963 ) performed a scale 
analysis which indicated that, in the absence of condensation, 
large-scale motions in the tropics are governed by the baro- 
tropic vorticity equation. Charney ( 1963 ) suggested that 
easterly waves arise from barotropic instability. The 
spectral analysis of Wallace and Chang ( 1969 ) indicate a 
conversion of kinetic energy from the mean flow to the 
disturbance . 
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A third possible energy source for these disturbances 
is latent heat of condensation. 

This mechanism is consistent with observations of warm 
core structures at upper levels (Chang et al. 1970). This 
heating has been parameterized by Ooyama (1964), Charney 
and Eliassen (1964), and Kuo ( 1965 ) in such a manner that 
the condensation heating is proportional to the convergence 
of moisture in the Ekman layer. This leads to a simple 
instability which is called Conditional Instability of the 
Second Kind (CISK). Yamasaki ( 1 9 6 9 ) studied the influence 
of CISK on linearized disturbances in an atmosphere with 
vertical shear. The solutions displayed many features that 
are observed in actual tropical wave disturbances. Holton 
(1971) has obtained similar results with a model that con- 
tains a moving heat source and vertical wind shear. These 
studies strongly suggest that the tropical disturbances are 
driven by condensation heating. It is known that the CISK 
mechanism must have a pre-existing disturbance to organize 
the convective elements before CISK can operate. This pre- 
existing disturbance could arise from barotropic instability 
as was demonstrated by Bates (1970). Kuo (1949) examined the 
stability characteristics of a barotropic zonal current and 
found that instability is possible only when the absolute 
vorticity gradient reverses sign in the region. 

The purpose of this thesis is to extend the work done 
by Williams et . al. (1971) to include cases with CISK heating, 
horizontal friction and vertical shear. Since the Williams 
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et al. (1971) paper is referenced many times, it will 
hereafter be referred to as WSN. 
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II. THE FORECAST EQUATIONS 



The quasi-geostrophic prediction equations will be used 
as the forecast equations following WSN. It is recognized 
that the quasi-geostrophic approximation is not very accurate 
near the Equator^ but Matsuno (1966) has shown with a baro- 
tropic model that quasi-geostrophic motions are identifiable 
near the Equator if the wavelength is not too long. 

Yamasaki ( 1 9 6 9 ) has found that the error is small at 20 
degrees latitude in a study with a baroclinic model. 

The quasi-geostrophic equations are simplified through 
the introduction of a two-level model which was also used 
by Bates (1970). The equations are linearized and solved 
numerically by the initial value technique. The growth rates 
and wave structures are then computed. 

A simple two-level model is constructed by dividing the 
atmosphere into four layers with a constant pressure differ- 
ential of Ap/2 (Fig. 1), numbered 0 to M from top to bottom. 
The boundary condition at p = p, is w = 0 , which prohibits 

Q 

vertical propagation of energy. Charney (1969) has shown 
that such energy propagation is unlikely. At p = p Q the 
boundary condition is w = w e where w e is the frictionally 
induced w at the top of the Ekman layer. Charney and 
Eliassen (19^9) derived the following expression for io e 



ur 



= -Pjjg (2J-1 sin 2 ct S| 



( 2 . 1 ) 
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where A g is the kinematic eddy viscosity, a is the surface 
inflow angle, and is the surface geostrophic vorticity. 



P 

p t 

P t + Ap 
p o 



w=0 





( 0=(0 

e 



Level 

0 

1 

2 

3 



Fig. 1 Two-level prediction model 



It has been shown by Holton et al. (1971) and Chang (1973) 
that this expression is reasonably valid as far equatorward 
as 10 degrees, but may be invalid when close to the Equator 
or south of the ITCZ. The surface geostrophic vorticity can 
be approximated by 

Cjj - = ^ ^3 * (2.2) 

The quasi-geostrophic vorticity equation in pressure 
coordinates with horizontal eddy diffusion is 

ft v2 * + fexvt-vcv 2 *) + e 0 || - f 0 If. - a h 7% . (2.3) 

Observe that when w = 0 and = 0, (2.3) becomes the baro- 
tropic vorticity equation which is sufficient to describe 
the main instability. 
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Apply this equation at levels 1 and 3 which gives 





(2.5) 



If it is desired to stop at some intermediate point, say the 
tropopause, rather than to include the entire atmosphere, 
smaller layers could be used. 

Next consider the quasi-geostrophic first lav; of thermo- 
dynamics in the form 



and T is the temperature obtained from the standard atmosphere. 

The heating follows the formulation of Ooyama (196*0, 
Charney and Eliassen (1964), and Kuo (1965), i.e., the 
condensation heating is proportional to the convergence of 
moisture in the Ekman layer. The form used was given by 
Ogura (1964) 





( 2 . 6 ) 



2f 

o 



where 




(2.6a) 
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i C p^ ae 

^2 ~ 2 . — ^ 

^ d e 3p 4 



(2.7) 



where n is a non-dimensional parameter. This relationship 
is used in both rising and subsiding air since a Fourier 
expansion in x is used; therefore the amplitude of heating 
is divided by two. 

When (2.6) is applied at level 2 the result is 

2 

a ’h+’f'o ApCpWp -u KQp 

^1-^3) + lkxV( x 2 - J ) • V(i{) 1 -4) 3 ) = — 2f — * (2,8) 



Define the following quantities 






4' 1 + 4» 3 
2 




which implies that' 

+ » 
* 



(2.9) 



( 2 . 10 ) 



( 2 . 11 ) 



( 2 . 12 ) 



Here ^ is proportional to the layer thickness and is there- 
fore a measure of the mean temperature in the layer. Use 
these definitions, add (2.4) and (2.5) and divide the result 
by two, which gives 
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It + lkxVi|; M .V(V 2 ^ M ) + lkxVijyV(V 2 <{> T ) 



+ 3 



d\j>. 
o 9x 



M 



f o“K 

2Ap 



V\ 



(2.13) 



For the second forecast equation, subtract (2.5) from (2.4) 
and eliminate with the use of (2.8) which gives 



9 ( V 2 -y 2 )if> T + lkxV^ M -V(V 2 -y 2 ) 1 4 T + IkxVijy V( V 2 ^) 



3t 



+ 3 



9 *T + f o w 4 



o 8x 2Ap 



-P 2 KQ 2 u 

2f + A H V j 
o 



where 



U 2 = 



2f 



P 



(2.14) 



(2 . l4a) 



These are the prediction equations for the model. To 

linearize the equations, the flow is separated into an east- 

west current that varies only in y, and a small departure 
from this flow. Thus it is possible to treat waves in x 
independently and the fields can be defined as follows 

= E(y) + A(y,t)cos kx + B(y,t)sin kx , (2.15) 

<4 T = F(y) + C(y,t)cos kx + D(y,t)sin kx , (2.16) 
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where k is the x wave number. The coefficients A through D 
are Fourier amplitudes of the disturbance and E and F are 
the coefficients of the mean zonal fields. 

Substitute the expressions for ^ and into (2.13) } 
separate the various sine and cosine terms, and neglect all 
products of the quantities A through D. Equate the coeffi- 
cients of the cosine terms which gives 






3 Bk - K(- — 7r - k 2 ) ( A - C) 
° 3y 




(2.17) 



The sine terms give 




k r 3E 9 2 A 3F 3 2 C 




+ — s- A + — 

3y J 3y 





(2.18) 
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Repeat the procedure for (2.14) for the cosine terms which 
gives 



§r ( i^ . Ck 2 - Cy 2 ) 
3t 9y 2 



D(k 2 + y 2 ) 

3y 3y^ 3y 3y^ 3y 



— BCk 2 - y 2 ) - ^|d - i^B] - B Dk 
3y 3y^ 3y^ ° 



+ (K - s) (■ 



3y' 



k^)(A-C) + A„(- — jj- 

3y 



2k 2 2% + A) 

3y 



(2.19) 



The results for the sine terms are 



3 (ifD . Dk 2 _ Dp2) . k C 9E 0(k 2 + ,2) + 3£ A(k 2 . p2) 

n by d by by 



1EA_^A + A fl + A c] + g oCk 



3y 3y 3y 3y 3y 



3y 



+ (K - S) (• 



3y‘ 



k^)(B-D) + Ajt(^— jjj- 

3y 



2 3 ^ 

2k^ i-4 + k D) , 

3y 



( 2 . 20 ) 



where 



fg a m ^ 
K = -2- (A 

2R0 o 



( 2 . 21 ) 
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and 



S 



g ri 9 0 / A Mx*s 1 

f | e | 3 P p ° V f o 



( 2 . 22 ) 



These are the same equations used by WSN except that in 
their equations F = 0, q = 0, and A^ = 0 . 

It is necessary to derive the computational equation for 
vertical motion beginning with 



03 



2 




KQ 2 

9t + U m 9x + V m 9y + 2f ^ * 

J o 



where 



u 



m 



1 

f G sy ’ 



v m = 



_ 1 _ 

f 3x 
o 



(2.23) 



(2.24) 



(2.25) 



Substitute equations (2.24), (2.25) and (2.16) into (2.23), 
collect terms, and again neglect all products of A through D 
which gives 



2 / , r 9C k ,9E n 9F n . 

u 2 Apa \ cos kx ^ 3t “ f Q ( 3y D “ 3y B ^ 



+ sinkx[f ^(ffc-ff A)] + f + ^- 



( 2 . 26 ) 
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Both the disturbance and the average vertical motion are 
given by this equation. The prediction equations and the 
vertical motion equations are now in the form to be used in 
the model. 
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III. BOUNDARY CONDITIONS AND FINITE 
DIFFERENCING SCHEME 

The finite difference scheme that was used is illustrated 
below with a sample variable : 



9M 

ay 



— (M 
2H u i+1 



M 



i-1 



), 



(3.1) 





2M ± + 



(3-2) 




2H 



3 [(M i + 2- 2 M i + r M i)- (M r 2 M i-i +M i-2 )] - 

(3.3) 



where i is the grid index and H is the distance between grid 
points . 

Centered time differences were used for all quantities 
except those involving friction and heating. The friction 
and heating terms were computed at time (t - At). Initially 
a forward time step was used. The second order equations 
for time tendencies were solved by the exact method of 
Richtmyer (1957, p. 101). 

Yanai and Nitta (1968a) studied the problem of using 
finite difference approximations to solve dynamic instabil- 
ity problems of non-divergent barotropic currents. They 
showed that for either symmetric or anti-symmetric zonal 
currents, the exact boundary conditions can be replaced by 
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placing rigid boundary conditions at a distance equal to 
one-half the width of the shearing wind belt. Therefore 
rigid boundaries are placed at y = -W/2 and y = W/2, where 
W is the total width of the computational field. Except 
for the E and F fields, all variables have the boundary 
condition A=B=C=D=0aty=± W/2. Second derivatives 
that are required on the boundary are set equal to zero. 
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Fig. 2 Mean wind profile U = -U tanh y/y Q 
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IV. WIND PROFILE AND INITIAL CONDITIONS 



WSN investigated the structure of tropical disturbances 
using a similar model and the wind profile (U = -U Q tanh y/y Q ). 
It is proposed to expand on that paper to include a sinusoidal 
profile for an initial disturbance and a thermal wind field 
which will be a mirror image of the mean wind field used 
by WSN. 

If the waves are to grow, the zonal wind profile must 
be barotropically unstable. Charney (1963) suggested that 
a shear zone could be formed by the transport of angular 
momentum from different regions. This would occur when the 
Intertropical Convergence Zone is shifted away from the 
Equator. The unstable profile chosen to represent this is 

Ui = U 3 = -U o tanh y/y 0 , (4.1) 

which is illustrated in Figure 2. To represent a thermal 
wind field the unstable profile is 

\J 1 = -U 3 = U Q tanh y/y Q . (4.2) 

Lipps (1970) has used the profile represented by (4.1) 
in an inviscid study of barotropic instability in the tropics. 
This wind field is anti-symmetric with respect to the 
Intertropical Convergence Zone, located at y = 0. 
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For the experiments which have the wind profile (4.1), 
a range of non-dimensional wave numbers (ky Q ) from 0.10 
to 3*00 will be used. This wide range of wave numbers is 
necessary because it is uncertain which wave number will 
be the most unstable mode. Those experiments with the (4.2) 
wind profile will be run over a range of 0.10 to 1.20. One 
experiment with the latter wind profile containing CISK 
heating was run over a wave number range of 0.10 to 3*00. 

tyjyj is introduced into the model as a sinusoidal distur- 
bance in y and ^ grows from zero as a result of the Ekman 
friction term. CISK heating is added which greatly increases 
the growth rates of all fields, particularly the short waves. 
Horizontal friction is also added. This represents a mixing 
of the smaller scale waves and smooths out perturbations 
in the vertical velocity field. Thus at time t = 0 



A = sin “Tj“ 
W 



( 4 . 3 ) 



B = 0 



( 4 . 4 ) 



C = 0 



( 4 . 5 ) 



D = 0 



( 4 . 6 ) 



E = f Q U q y Q ln[cosh (y/y Q )] or 0 (4.7) 



F = -f U y ln[cosh (y/y )] or 0 (4.8) 

o o o o 
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V. COMPUTATIONAL PROCEDURE 



The following constants were used in all computations 
except those related to heating and friction which were 
zero if the heating and/or the friction are neglected: 



p o 


= 


100 centibars, 


y o 


= 


400 or 200 or 100 kilometers, 


W 


= 


4000 kilometers. 


U 

o 


= 


10 meters per second, 


f 

o 


= 


_5 

5 x 10 per second, 


c 


= 


0.4 meters per second squared per 
centibar squared, 


B 


= 


0.0 or 2.29 x 10”'*''*' per meter per second 


Ap 


= 


50 centibars, 


A 

e 


= 


10 meters per second squared, 


a h 


= 


5 

1 x 10 meters per second squared. 


a 


= 


22.5 degrees, 


2 

y 


= 


—12 

2.49 x 10 per meter squared, 


H 


= 


50 or 25 or 12.5 kilometers. 


At 


= 


1.0 or 0.5 or 0.25 hours. 


n 


= 


2.0 or 1.0, 


36/3p 


= 


-0.6 degrees per centibar. 



A forecast period of 142 days was used so that the most un- 
stable mode would have reached the maximum growth rate. For 
most cases, it was found that the growth rate at 110 days 
was the same as the growth rate at the longer period. 
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WSN introduced a random disturbance into the ^ field. 

It was found that by using a sinusoidal profile for the 

disturbance, similar growth rates were obtained except in 

the cases where heating was added. The sinusoidal profile 

was used as the disturbance field for all experiments with 

variations of U, 3, CISK heating, horizontal and surface 

friction, and the non-dimensional wave number (ky ). Growth 

o 

rates were then computed by assuming the amplitude could 
be written as 



a i = a o e xp( nU 0 /y 0 ) t l» 



(5.1) 



a 2 = a 0 ex p( nU 0 /y Q ) 



(5.2) 



where n is the growth rate and a Q is a constant. If the 
ratio a 2 /a-^ is formed, the growth rate can be shown to be 



n 



y Q ln(a 2 /a 1 ) 

U o (t 2 - V 



The dimensional doubling time is 



(5-3) 



T 



y 



In 2 
o 

nU o 



(5.*0 
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VI . RESULTS 



A series of experiments were conducted over a wide 
variety of conditions. The total distance over the grid 
was JJ000 km. A grid mesh size of 50 km and a time step of 
1.0 hours was used for all experiments except two cases. 

In one case a grid mesh of 25 km with a 0.5 hour time step 
was used and in the other case a grid mesh of 12.5 km and 
a 0.25 hour time step was used. 

The more important experiments are given in Table 1 as 
a function of U, 3, heat, friction, and ky o - The growth 
rate n is also given. The table shows whether or not CISK 
heating and horizontal friction were included for each 
experiment. All experiments given in Table 1 included 
surface friction. For each experiment, calculations were 
made to determine the non-dimensional growth rate as a 
function of the non-dimensional wave number ky so that the 
results could be compared to those obtained by WSN. The 
value of ky Q that was associated with the maximum growth 
rate is shown in the table. 

It should be mentioned that the value used for the non- 
dimensional heating parameter n was too large. A more 
realistic value for n would have been 1.0 instead of 2.0. 
Test runs were made for experiments 1 through ^ with 
n = 1.0 and in all cases, the growth rates were reduced to 
about one-third the value obtained with the larger value 
of n. 
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Table 1. Compilation of numerical experiments 



Exp 


u 


e 


n 


Heating 


Friction 




1 


°1 ■ 


0 


0.770 


Yes 


No 


0.40 


2 


“l * “3 


2. 29x10“ 1:L 


0.730 


Yes 


No 


0. 40 


3 


U 1 = U 3 


0 


0.520 


Yes 


Yes 


0.40 


4 


“l = «3 


2.29xl0 -11 


0.520 


Yes 


Yes 


0. 40 


5 


U 1 ’ ' U 3 


0 


0.283 


No 


No 


0.65 


6 


U 1 ’ - U 3 


0 


0.790 


Yes 


No 


1.10 



Experiment 1 was run with the (4.1) wind profile and 
included CISK heating and surface friction. The beta term 
and horizontal friction were excluded. Garcia (1956) found 
that the shear zone profile (4.1) is barotropically un- 
stable when 0 < ky Q < 1. Betchov and Criminale (1967) com- 
puted the eigensolutions for this profile throughout this 
range and found the maximum instability at ky Q = 0.45 which 
agreed with the results of WSN. The wave number of maximum 
instability ( 0 . 4 0 ) agrees quite closely with that value but 
the growth rate (n = 0.770) is much greater than the value 
obtained by WSN (n = 0.153) even when the smaller value of 
h was used. Figure 3 contains the growth rate profile for 
the full range of wave numbers. The maximum growth rate 
corresponds to a doubling time of 0.42 days and a wavelength 
of 6280 km. The sharp irregularities between 0.10 and 1.00 
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arise because the wave number increment in that region is 
0.10 as compared to 0.50 between 1.00 and 3.00. In the 
non-heating cases of WSN, there was no growth for ky Q > 1, 
but in experiment 1 there is large growth which must be due 
in part to the CISK process. The Figure shows that the 
growth is nearly independent of the scale. Figure 4 contains 
the amplitudes of ^ and oo 2 for this experiment. ^ and ip^ 
are not shown for clarity purposes as they have the same 
magnitude as ip^ is slightly larger and ip^ slightly 

less than ip^ but the shape is identical. The units are 
arbitrary since the basic equations are linear. However, 
u 2 is not plotted on the same scale as ip . All quantities 
show double maxima with peaks at y = ±1.2. These results 
are somewhat similar to those obtained by WSN without 
heating but their results showed the minimum between peaks 
as only slightly less than the maximum whereas in Figure 4, 
the minimum is less than one-half of the maximum. Also the 
14^, ^3 j and i4ip fields are very close together but WSN 
showed 14^ to be predominate. The closeness of the ip fields 
roughly implies that the vertical mean streamfunction is 
nearly zero so the fields are basically of a thermal charac- 
ter. Figure 5 contains the phase angles for 14^, '4 T > 

and oo 2 for this experiment as functions of y. ^ and ip^ 
are very nearly coincident at the outer edges and 14 = ip ^ 
between y = ±2 although only 14^ is shown in Figure 5. ^3 

is approximately 180 degrees out-of-phase from 14^ and i4rp 
which indicates the pressure trough at the surface is replaced 
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Pig. 3 Growth rates for experiment 1. 
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Fig. 4 Amplitudes of i4 T , and w 2 for experiment 1. 
tfj,p and ta ^ are on different scales. 
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Pig. 5 Phases’ of if» T , and w 2 for 

experiment 1 in degrees. Each dependent 
variable is expressed in the form 
Pcos(kx-S), where P is the amplitude 
and 6 is the phase angle. 
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by a ridge at the upper level. This is consistent because 
the Ekman effect and heating are tied to the field at level 
3 which produces heating above surface low pressure systems. 
Unlike the results found by WSN for the adiabatic case. 

Figure 5 shows the horizontal tilt of the disturbance to 
be a minimum at y = 0 and a maximum at y = ±(1 to 4). The 
general sense of the tilt is northeast but larger than that 
observed in easterly waves. Also the tilt of the o> 2 field 
fluctuates considerably. The small perturbations, particu- 
larly in the u> 2 field, suggest that there is not a clearly 
defined unstable mode in y. This does not agree with the 
results of WSN which had no heating. It is interesting 
to note that the fields are very smooth but have large 
gradients except for w 2 which shows a radical behavior. 

Experiment 2 is the same as experiment 1 except the 
beta term is included. The results in general are the same. 
Figure 6 shows the growth rate profile. The magnitude of 
the growth rates is almost the same although a small decrease 
in the maximum is noted (from 0.770 to 0.730) and the 
perturbations betwen ky Q = 0.10 and 1.00 have been somewhat 
smoothed. The maximum growth rate corresponds to a doubling 
time of 0.44 days and a wave length of 6280 km. Figure 7 
contains the amplitudes for experiment 2 as a function of 
y. As in experiment 1, only ^ and co 2 are shown with <4rp 
and co^ on different scales. Again there is a double maxima 
at y = ±1.2 but the profile is asymmetric with the larger 
maximum in the shear zone to the north. This asymmetry is 
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Fig. 6 Growth rates for experiment 2. 
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pig. 7 Amplitudes of i^.and Wj for experiment 2. 
4 > t and Wg are on different scales. 



33 




Fig. 8 Phases of tyy ^ T > and u 2 for 

experiment 2 in degrees. Each dependent 
variable is expressed in the form 
Pcos(kx-6), where P is the amplitude 
and 5 is the phase angle. 



3 



expected since the absolute vorticity gradient is not an 
odd function of y. VJSN found that the shift of symmetry was 
to the south while Figure 7 shows a shift to the north. 

The phase relationships of 4^> and w 2 shown in 

Figure 8 are very similar to those in experiment 1 and 
again only 4^ is shown between y = ±2. 

Experiment 3 was run with the (4.1) wind profile and 
included CISK heating, surface friction, and horizontal 
friction. The beta term was excluded. The horizontal 
friction was added in an attempt to smooth out the vertical 
velocity (u^) pattern. Figure 9 shows the growth rate 
profile for experiment 3. The maximum growth rate has been 
reduced considerably (from 0.770 to 0.520) and the growth 
falls off at ky Q >_ 2.00. The maximum growth rate corresponds 
to a doubling time of 0.62 days and a wavelength of 6280 km. 
Note that there are no perturbations in the profile. Figure 
10 contains the amplitudes of 4'rp and for this experiment. 
As in experiment 1, only 4^ and w 2 are shown and are on 
different scales. '4^ is slightly larger and 4^ slightly 
smaller than 4^ • It is interesting to note that the largest 
amplitudes are located near the boundaries for all fields 
whereas in the case without friction (experiment 1), the 
maximum amplitudes are in the central regions. Perhaps the 
large scale eddies favored by the horizontal friction are 

v 

heavily damped by the horizontal shear between y = ±1. 

Figure 11 contains the phases of 4^> 4^ and a) 2‘ ^T no ^ 
shown but it is everywhere equal to 4^* Again ij'ip and ip^ 
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Fig. 10 Amplitudes of and w 2 for experiment 3- 
4>rp and w 2 are on different scales. 
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Pig. 11 Phases of 4^, and oo 2 for experiment 

3 in degrees. Each dependent variable 
is expressed in the form Pcos(kx-5), 
where P is the amplitude and 6 is the 
phase angle. 
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are out-of-phase but all fields are smoother than in experi- 
ment 1. Where there are peaks in the amplitude fields, the 
phase angles of and ^ are the same which tends to 
suggest that warm air is rising and cold air is sinking. 

Again there is an eastward tilt in the <4 fields. The io 2 
field tilt is opposite to the fields but it does correlate 
with the iK field where the amplitude is a maximum. It 
appears the scale in y is about 2y Q . 

Experiment 4 is the same as experiment 3 but the beta 
term is included. The beta term has very little effect 
except to make the solutions asymmetric. Figure 12 shows 
the growth rate profile which is the same as Figure 9 
except for the slight perturbation between the wave numbers 
0.40 and 0.90. The maximum growth rate corresponds to a 
doubling time of 0.62 days and a wavelength of 6280 km. 

Figure 13 contains the amplitudes of i^rp and co 2 • As in 
experiment 3, ^ and are not shown. The only change noted 
is a slight increase in the maximum to the north and a 
decrease in the maximum to the south. Figure 14 contains 
the phases of and w 2 and again i4fp is everywhere equal 

to The w 2 and ^ correlation is the same as in experiment 

3. 

Experiment 5 was run with (4.2) wind profile and surface 
friction. Heating and horizontal friction were excluded. 

This profile is a mirror reflection of the profile, that 
is -U^. This profile is similar in general to the wind 
profile in the Intertropical Convergence Zone where there 



39 



is cyclonic shear in the lower region and anticyclonic 
shear in the upper region. Since the ITCZ is driven by 
heating and convergence, this is the type of mean wind 
pattern that is expected. Three cases were run although only 
one is shown in the table. Case 1, which is summarized in 
Table 1, had a grid mesh of 50 km, a y Q value of 400 km and 
a time step of 1.0 hours. Case 2 used a grid mesh of 25 km, 
a y Q value of 200 km and a time step of 0.50 hours. Case 3 
used a grid mesh of 12.5 km, a y Q value of 100 km and a 
time step of 0.25 hours. Case 1 was run over the range of 
wave numbers from 0.10 to 1.20 and the maximum unstable 
mode was found at ky Q = 0.65* Cases 2 and 3 were run only 
at wave number 0.65. Figure 15 shows the growth rate profile 
for experiment 5, case 1. The maximum growth rate at ky Q = 0.65 
was 0.283. The doubling time was 1.13 days with a wavelength 
of 3465 km. For case 2, the growth rate obtained was 0.188 
with a doubling time of 0.85 days and a 1930 km wavelength. 

For case 3 a. growth rate of 0.162 was found with a doubling 
time of 0.50 days and a 967 km wavelength. For cases 2 and 
3, ky Q = 0.65 may not be the maximum growth rate. This shows 
that as the scale (y Q ) is decreased, the baroclinic process 
is suppressed. This is consistent with baroclinic instab— 
bility theory. It is noted that in case 3 , the growth rate 
is the same as that obtained by WSN in their experiment 1. 
Figure 16 contains the amplitudes of 4^, and w 2 for case 
1. For clarity, 4 ^ i s n °t shown but it is slightly less 
than i|/ . These amplitudes are very similar to those of WSN 
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Fig. 13' Amplitudes of and for experiment 
4. ^ and are on different scales. 
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Pig. 14 Phases of i|^> and w 2 for ex P erlment 

4 in degrees. Each dependent variable 
is expressed in the form Pcos(kx-S), 
where P is the amplitude and 6 is the 
' phase angle . 
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in their Figure 4A. The amplitudes for cases 2 and 3 are 
not shown but they are much the same as Figure 16 . Figure 
17 contains the phases of i|^, ij^, and as functions 

of y. and are approximately 165 degrees out-of-phase. 
This implies a baroclinic conversion of kinetic energy. 
However, the tilt is opposite to the shear which indicates 
a barotropic conversion of kinetic energy. The tilt is 
westward with height north of the ITCZ and eastward with 
height south of the ITCZ which is consistent with rising 
motion and unsettled weather that is often observed east 
of the surface trough axis in easterly waves north of the 
ITCZ. These results are also consistent with those of 
WSN. 

Experiment 6 was run with the (4.2) wind profile, CISK 
heating and surface friction. The beta term and horizontal 
friction were excluded. Figure 18 shows the growth rates 
for this experiment with a maximum unstable mode at ky Q = 1.10 
corresponding to a doubling time of 0.4l days and a wave- 
length of 2285 km. The growth rate of 0.790 is similar to 
that of experiment 1 and again the large value was due to 
the large value of n. The perturbations between wave numbers 
0.10 and 1.10 are even larger than in experiment 1 indicating 
that there is no clearly defined unstable mode in y. The 
amplitudes have a pattern similar to Figure 4. The phases 
for this experiment tilted northeast as in the other cases. 



44 




Pig. 15 Growth rates for experiment 5, case 1. 
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Fig. 16 Amplitudes of i|; , and w 2 for 

experiment 5, case 1. The quantities 
and ipy have the same scale; is 

on a different scale. 
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Pig. 17 Phases of ^ T > and w 2 for 

experiment 5, case 1 in degrees. Each 
dependent variable is expressed in the 
form Pcos(kx-6), where P is the 
amplitude and 6 is the phase angle. 
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VII. CONCLUSIONS 



In this thesis, numerically computed barotropically un- 
stable disturbances are compared with tropospheric waves 
that occur in the tropical easterlies. A two-level quasi- 
geostrophic model is used with an attached Ekman layer. The 
most unstable solutions are obtained numerically by the 
initial value procedure. Two basic shear zone profiles are 

used: (1) U, = U = -U tanh y/y and (2) U n = -U_ = U tanh y/y . 

1 3 o JJ o 1 3 o J J o 

Calculations were made in which CISK heating and horizontal 
friction were added. 

Certain features of the predicted wave structures do 
correspond to observed tropical wave disturbances. The tilt 
of the disturbance in a sense opposite to the horizontal 
wind shear that was observed by Wallace and Chang (1969) is 
present in all the experiments. All experiments which contain 
heating have the warm structure which was observed by Chang 
et' al. (1970). The presence of rising motion to the east 
of the wave axis which was noted by Yanai and Nitta (1967) 
is observed in experiment 5. 

The growth rates obtained with heating are unrealistically 
large, but it was shown that these could be greatly reduced 
by choosing a smaller value of the parameter p . It was also 
noted that when heating is present there is no truely 
preferred scale of instability. This has been observed by 
Chang (1971). It was pointed out in the Introduction that, 
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while mature tropical wave disturbances may be driven by 
CISK heating, they may originate from barotropic instability 
(Bates 1970). If so, then the disturbance should retain 
certain features (latitude structure, wavelength, etc.) 
of the barotropic ally unstable modes. This would provide 
a selection process which does not occur if CISK heating 
is present from the start. 

The studies of this thesis should be continued but a 
more complicated model is needed. Proper parameterization 
of the heating term is necessary. A primitive equation 
model which includes a better treatment of the boundary 
layer should be used. 
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